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Quantum vortices in the multi-component Bose-Einstein condensation (BEC) are investigated 
theoretically. It is found that three kinds of the vortex configurations are possible and their 
physical properties are discussed in details, including the density distribution and the spin tex- 
ture. By using the Bogoliubov theory extended to the three component BEC, the collective 
modes for these vortices are evaluated. The local vortex stability for these vortices are examined 
in light of the existence of the negative eigenvalue, yielding a narrow magnetization window for 
the local intrinsic stable region where the multi-components work together to stabilize a vortex 
in a self-organized way. 

KEYWORDS: Vortex, Three component Bose-Einstein condensation, Gross-Pitaevskii equation, Bogoliubov equa- 
tion, Collective modes 
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§1. Introduction 

There have been much attention focused on Bose- 
Einstein condfmsation (BEC) realized in alkali-metal 
atom EasjesJd' h\3 In particular, since two experimental 
groupsa'BBl3'El> have succeeded in producing quantum 
vortices in a rotating trapped BEC in 87 Rb system, many 
researchers investigate various aspects of this quantized 
vortex in the hope to obtain the better understanding of 
physical nature of superfluidity in this dilute weakly in- 
teracting Bose system. Here we have a better controlled 
mathematical tool, namely, Bogoliubov theory to treat 
it, compared with superfluid 4 He where no established 
microscopic theory to treat this strongly interacting Bose 
condensed system is knowji 

In a series of our papersollJ : E]E3 ) we have been study- 
ing the problem of the intrinsic instability of a vortex in 
BEC. We point out the existence of the negative eigen- 
value in the Bogoliubov excitation spectrum^ demon- 
strating that a quantized vortex in BEC confined by a 
rigid wall is intrinsicalbz_.unstable and corresponding to 
the Rokhsar instability^ due to the spiraling-out of a 
vortex core in a harmonically trapped BEC. Indepen- 
dently, the existence of the negative value is pointed by 
Dodd et alLji 1 In order to arrest this instability and to 
stabilize it, we propose several practical devices feasible 
experimentally, such as sending a laser light at the core 
center to raise the potential which moves up the whole 
Bogoliubov spectrum as a result, or raising temperature 
and increasing the non-condensate fraction relative to 
the condensation which acts as an extra potential^ A 
similar recovering method is proposed by Doill^ who con- 
sider the two compopent BEC. We also examine the ex- 
ternal rotation effectE2f on the vortex stability and deter- 
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mine_the several characteristic critical rotation frequen- 
cies,^ including the lowest critical frequency which sig- 
nals the first entry of a vortex unon increasing the forced 
rotation. This particular valueEP is by far nearer to the 
observed stirring frequency in-nthe recent vortex creation 
experiment by optical spoonHLrB) than that estimated 
by the conventional Thomas-Fermi approximation. This 
coincidence further encourages us to perform the present 
study along this line. 

So far in all the theories and experiments mentioned 
above the internal degrees of freedom in an atom is frozen 
by the external magnetic field serving as a confining po- 
tential. By means of optical confinement, all the hy- 
perfine states with F=l in 23 Na active and three com- 
ponents with tof = 1,0,-1 are simultaneously Bose- 
condensedjij) giving rise to a new state of mattex: a BEC 
characterized _hy a vectorial order parameter .EIP Ohmi 
and Machida,E3 and HoES* independently introduce the 
basic Hamiltonian for describing this vector BEC by ex- 
tending the Bogoliubov framework and study the funda- 
mental properties of this interesting vector BEC, point- 
ing out the richness of the topological defect structures 
such as the ^-vector textures, domain wall structure, etc. 
These are very analogous to the 3 He superfluid A and B 
phases.EP Again as the basic interaction in these super- 
fluid Fermions is strong, we have to resort to the strong 
coupling theory to describe it. The existing theories for 
superfluid 3 He is quite phenomenological and difficult to 
derive microscopically. In contrast, our theory here is 
quite microscopic and there is no room to introduce phe- 
nomenological parameters in it. 

Guided by this analogy, we are going to_ study here 
the vortex structures in more details. YipElP gave the 
first study on a vortex in the vectorial BEC, classifying 
the possible vortex structures as functions of the interac- 
tion strength of the spin channel and the magnetization 
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of the system. Unfortunately, some of the vortices in- 
troduced by Yip (I and IV phase in ref. ^l|) are simply 
not allowed because the fundamental phase requirement 
is violated for these particular vortices (See eq. ([To|)). 
Lconhardt and VolovikESl-ciropose the half-quaptum vor- 
tex (Alice string). StooS and Marzlin et a£-$ exam- 
ine independently the skyrmion creation in the vectorial 
BEC .-JiuaSr. here . continuing our study on the vectorial 
BECjla'Ea'E^'EZl'ES) we perform microscopic calculations 
of the possible vortices in harmonically trapped vector 
BEC and examine the intrinsic stability problem in con- 
nection with our intrinsic instability of a vortex in the 
scalar BEC. 

After giving the Gross-Pitaevskii equation for the vec- 
torial BEC with the three component order parameters 
in next Section, we examine the fundamental phase re- 
quirements for the order parameter components when 
they form a vortex. We enumerate several possible vor- 
tices in a cylindrical symmetric case in §3. In §4 each of 
these possible vortices is studied in details: their physical 
properties such as the density distributions of the three 
components and spin configuration around a vortex core. 
Evaluating the excitation spectrum for each vortex form 
by solving the Bogoliubov equation extended to the three 
components, we examine the stability condition of each 
vortex as a function of the total magnetization in §5. 
The last Section is devoted to summary and discussions. 

§2. Gross-Pitaevskii Equation for Three Com- 
ponents 

We treat a system of the Bose condensates with the 
internal degrees of freedom F = 1 so that the order pa- 
rameter is characterized by the spins: mp — 1,0,-1. 
We start with a general Hamiltonian which fully consid- 
ers system's symmetries in real space and spin space: 



cedure for a BEC system without spin, the Gross- 
Pitaevskii equation extended to the three component or- 
der parameters is obtained as 

[jfr(r)-/i !: + ffn ^|^| 2 j<% 
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These coupled equations for the j-th condensate <pj (j = 
0, ±1) are used to calculate various properties of a vortex 
in the following. 

§3. Phase Requirements of Condensates 

The two-dimensional disk system is assumed and the 
coordinate system r = (r, 9) is used. The wavefunctions 
are divided into the radial part <p and phase part 7 as 



7j -(0)=exp[i(a j +&0)] 



(6) 
(7) 



where j = 0,±1. The wavefunctions </>j (r) are real. Tem- 
porarily we assume 4>j(r) > 0. 

The total energy of the system is given by 
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The phase coefficients jj of the condensate wave function 
are determined such that the energy of the g s term in 
integrand of eq. (8) is minimized. This is because the 
other terms in eq. (8) are not affected by the choice of 
phases. The condition to minimize the g s term: 
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leads to 



m 3 

The subscripts are a = (x,y,z) and i,j,k,l = (0, ±1) 
corresponding to the above three species. The chemical 
potentials for the three components; /ij (i — 0,±1) sat- 
isfy fii = 2/iQ. The scalar field V(r) is the external 
confinement potential such as an_optical potential. 
Following our previous paper,[3 and the usual pro- 
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±1. 
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where the upper (lower) sign is used when g s < (g s > 
0). This condition is rewritten in terms of a and (3 as 

2ao = oti + OL—i + nn (11) 

2/3 = + /3_ x (12) 

where n is an integer. We take ao = in the following. 
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Thus the phase factors of the wavefunctions are generally 
written as 



7i 
7o 
7-i 



exp[i(a + (39)] 

1 | exp[i(3 9] (13) 

±exp[-i(a + (39)} 



with 



and 



P = p 1 -po=Pa-P- 
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In the followings we exclude the vortex states with 
winding number /3j larger than 1 simply because the 
smaller winding number, the more stable. The possible 
combinations for (3i are (/3i,/3o, /3-i) = (1, 0, —1), (1, 1, 1) 
and (1, 5,0). Because the value 5" does not allowed in 
this cylindrically symmetric system, this means that this 
component vanishes. We write (1, |,0) as (l,none, 0) 
in the following. We call this vortex state the Alice 
state. JIhe Alice string discussed by Leonhardt and 
VolovikE3 has this combination of the winding numbers 
for the three components, although the Alice string has 
another condition for the amplitudes. The trivial com- 
binations which are equivalent to the above ones, i.e. 
(-1,0,1), (-1,-1,-1), (-l,none,0) and (0,none,±l) 
are excluded. 

The sign of 7_i in eq. ( [13] ) is determined by the sign of 
g s as shown in derivation of eq. ([l3|) . When g s < (fer- 
romagnetic case), 7_i = + exp[— i(a+(39)]. When g s > 
(antiferromagnetic case), 7_i = — exp[— i(a + (39)]. We 
redefine <p-i and 7_i for simplicity of our calculation: 



-i( r ) > (ferromagnetic case) 
-i( r ) < (antiferromagnetic case) 



(16) 



7_i =exp[-i(a + /30)] (17) 

The Gross-Pitaevskii equation eq. (5) for the real wave- 
functions 4>j (r) becomes 

-ft 2 f d 2 Id ((3 +j(3) 2 
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The phase a does not appear here. This means that 
4>j(r) are independent of a. 

We have done the extensive calculations to solve the 
above equations. To concentrate on a long straight vor- 
tex, we assume that the system is uniform along the z 
axis as mentioned. Since in the ferromagnetic case the 
phase separation occurs, the results for a disk shape here 
may be appropriate to this case. 

The properties of F = 1 condensate is classified into 
ferromagnetic and antiferomagnetic depending on sign 



of g s . To represent antiferromagnetic (g s < 0) and 
ferromagnetic^ < 0) case, g s — +0.1 and g s — —0.1 
are used as examples. 

The actual numerical computations are performed un- 
der the conditions: The mass and the scattering length 
of atoms are to — 3.81 x 10~ 26 kg and clq ~ 2.75nm. 
The area density of particle is 2 x 10 4 (/im) _1 . These 
are appropriate for those of Na atom. The other scat- 
tering length a-i is defined so that <7 S becomes ±0.1g n . 
The trapping frequency is uj/(2tt) = 200Hz. The energy 
is scaled by the trap unit Klo. The confining potential is 
given by V(r) = ^mui 2 r 2 . 

We calculate a set of the physical quantities under 
given parameters such as g s , the relative polarization 
M/N (see below for the definition) which is determined 
by adjusting the chemical potentials, and the combina- 
tion of winding number /3's. We vary the relative polar- 
ization from to 1. Using symmetries, this covers all the 
possible value of M/N. Both ferromagnetic (g s = —0.1) 
and antiferromagnetic (g s = 0.1) cases are comparatively 
studied in the following. 

§4. Three Kinds of Vortices 

We calculate the three kinds of the vortex states for 
various magnetizations. The number of atoms and the 
magnetization are given by N — J7=i — 1 an< ^ ^ ~ 
J2i=i -1 where Ni — J" |t^i(r)| 2 dr are the number 
of atoms in the spin state i. The range of the relative 
polarization is -1 < M/N < 1. As for the (1,0,-1) 
and (1,1,1) configurations, it is sufficient to consider the 
range < M/N < 1. 

It is convenient to describe the condensates in terms of 
the three wavefunctions 4> a (a — x, y, z) where the spin 
quantization axis is taken along the a direction: 








(19) 



Let define the to- and n-vectors as 

m = (m x , m y , m z ) = Re((f> x , <f> y , 4> z ) (20) 

n = {n x ,n y ,n z ) = lm((f> x , <j) y , <j> z ). (21) 

where m and n are real vectors. The Z-vector which 
points the direction of the local magnetization is 1 = 
m x n. The corresponding unit vector is 1. 

4-.1 Alice states 

When the phase configuration of the condensates is 
the Alice type: 




(22) 



the Z-vector is always parallel to the z-axis: l x = l y = 0, 
l z = 4>\ — ((pLi- The zero component of the condensate 
4>o is always zero and the system is the 2-component, 
consisting of 4>i and 4>—\. 
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When the interaction is ferromagnetic, it is conceivable 
that one of the two components vanishes and the con- 
densate reduces to the single component system. The 
system hardly has intermediate magnetization. How- 
ever, when the +1 component carries the phase wind- 
ing, the other component may appear at the center of 
the vortex. In fact, as shown in Fig. 1(a) where the 
density profiles of the Alice state are displayed, the 
core region of the +1 component is filled by the re- 
maining component -1. The range of magnetization is 
M/N = -1, and from + 0.967 to + 1. In the inter- 
mediate range —1 < M/N < +0.967, the 0o compo- 
nent, which is not allowed in this cylindrically symmetric 
system, is likely to appear due to spin interaction term 
eq. (9) and cylindrical symmetry solution of GP equation 
does not obtained.Er 

When the interaction is antiferromagnetic, the ampli- 
tude ratio of the two components may vary continuously, 
depending on the magnetization. Figure 1(b) shows the 
case when M/N = 0.531. For all ranges of M/N, the 
density of the +1 component has a peak at r ~ 3fxm. 



4.2 (1, 0, -1) vortex 

In this vortex, the condensate is described by 

(f>t \ / Mr)e iW ° +a) 
0o = <f>o(r) 
6-i / V 0_ 1 (r)e- i (' 3e+Q ) 



(23) 



The winding number of 4>o is zero and the component is 
only the component of the condensates which can exist at 
the vortex center. The magnetization of the condensate 
can vary continuously in this configuration of the winding 
number. The Z-vector is given by 

6i(r) + <f>-i(r) 
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— sin (a 



p6) 
f P0) 



024) 



(25) 



Figure 2 shows the l x and l y configurations in the (x, y) 
plane for the two kinds of the disgyration with cither 
= 1 or /3 = -1. 

Figures 3(a) and (b) show the particle density profiles 
along the radial direction in the ferromagnetic case. It 
is seen from Fig. 3(a) that the (f>o component fills in the 
core region, thus there is no appreciable particle number 
depleted region near the center. Notice that </>_i(r) < 
when the interaction is ferromagnetic. 

The corresponding antiferromagnetic case for the 
(1,0,-1) vortex is displayed in Fig. 4(a) and (b). Here it is 
also true that the large number of the c/)q component oc- 
cupies the core region, pushing out 4>i and <j>-\. The total 
particle number almost smoothly decreases outwardly. 




Fig. 1. The density profiles of the condensates are displayed 
along the radial direction and the /-vector in the Alice type vor- 
tex is shown for the quadrant in the (x, y) plane. The vertical 
axis is the density of the condensates. The upper bold curve 
shows the total density n(r) and the thin curves show the den- 
sity of each component n;(r). The arrows on the {x,y) plane 
are the Z-vectors. The small circles are perpendicular to each 
Z-vector a nd the sh ort rods show m-vectors. They are depicted 
at r = \J x 2 + y 2 = 0,1,2,3,5, and 7 (im. (a) Ferromagnetic 
case. The magnetization M/N = 0.967. (b) Antiferromagnetic 
case. M/N = 0.531. Note that the intermediate magnetization 
like this is not possible for the ferromagnetic case. 
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Fig. 2. The bold arrows show x and y component of Z-vectors 
when winding number combination is (1,0,-1) for the ferro- 
magnetic case (a = 0). The Z-vectors show cross disgyration 
when fi = 1 and radial disgyration when /3 = — 1. For the an- 
tiferromagnetic case the direction of the Z-vectors in the (x,y) 
plane is reversed. 
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Fig. 3. The density profiles of the condensates and the i-vector 
when the interaction is ferromagnetic and the winding number 
combination is (1,0, —1)- The vertical axis is the density of the 
condensate. The upper bold curve shows the total density n(r) 
and the thin curves show the density of each component n;(r). 
(a) M/N = 0.612. (b) M/N = 0.994. 



4-3 (1,1,1) vortex 

All of components have the winding number 1 in this 
vortex, namely, the wavefunction is written in the form: 





(26) 



It is shown that the /-vector does not depend on 9. 

^o^HcosM \ 
]= I foil^^-a) . (27) 

For the ferromagnetic case, this reduces to single com- 
ponent BEC. Figure 5 shows the typical results for this 
case, namely, as all the components have non-vanishing 
phase winding, the densities vanish at the vortex core 
center. The total density distribution looks like a vortex 
in the scalar BEC. For the antiferromagnetic case, the 
/-vector given by eq. (|27|) vanishes when \<p-i\ = \<f>i\. 



§5. Local Stability of Vortex State 

5. 1 Formulation 

In order to determine the intrinsic stability condition 
for each vortex state without the external rotation, we 




Fig. 4. The density profiles of the condensates and the i-vector 
when the interaction is antiferromagnetic and the winding num- 
ber combination is (1, 0, —1). The vertical axis is the density of 
the condensate. The upper bold curve shows the total density 
n(r) and the thin curves show the density of each component 
n;(r). (a) M/N = 0.474. (b) M/N = 0.952. 




Fig. 5. The density profile of the condensate and the i-vector for 
the ferromagnetic case of the (1, 1, 1) vortex. The magnetization 
M/N = 0. The vertical axis is the density of the condensate. 
The upper bold curve shows the total density n(r) and the thin 
curves show the density of each component rii(r). 
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give the formulation for evaluating the excitation spec- 
trum of the three component BEC, which is based on the 
Bogoliubov theory extended to take into account the spin 
degrees of freedom. The collective excitation spectrum of 
the condensate, whose static properties are determined 
by the Gross-Pitaevskii equation eq. (5), is a solution of 
the Bogoliubov equations given by 



3 



{AijUqirJ) - BijV q {v, j)} = e q u q (r,i), (28) 



£{B£«,(r, j) - A* jVq (r,j)} = e q v q (r,i) (29) 



where 



Aij = In r id,, - fiiSij 

+g n I ^2 l^fel 2 ^. 



'9az2z2[( F a)ij(F a )kl4>k ( l ) l 
a hi 

+(F a ) u (F a ) kj (ttfr)}, 



(30) 



Bij = g n <i>i<t>3 + 9s ^2 ^2 [{F<x)ik(t>k{F<x)ji(t>i] , (31) 

a kl 

u q (r 7 i) and v q (r,i) are the <?-th eigenfunctions with the 
spin i and e q corresponds to the g-th eigenvalue. 

We are treating a cylindrical system r = (r, 8, z). The 
relevant physical quantities are written in the following 
form up to the phase factor: 

u,(r, j) = u q (r, j) e Kie+0o)e e i(a+pe)3\ (32 ) 

Vq (r, j) = v q {r,j)e l{q »- Me e- l{a+l5e) \ (33) 

with the angular momentum qg. (The momentum q z 
along the z direction is ignored.) 

5.2 Stable range 

When all of the excitation levels bffiS the positive 
eigenvalue, the vortex state is stable Jj'EHP namely, the 
vortex is locally stable in the energy configuration space. 
Thus the stable region of the magnetization for a given 
vortex is determined by evaluating the whole excitation 
spectrum of the above Bogoliubov equations: eqs. (28) 
and (29). We calculate them to check if there exists 
the negative eigenvalue. This procedure yields the stable 
magnetization region. It should be noted that the lowest 
level at the particular angular momentum qg = —1 in 
the excitation spectrum often becomes negative, signal- 
ing the local instability-oLa vortex as discussed in the 
single component BEC.BlEJ The positive spectrum does 
not mean that this vortex is stable globally relative to 
other states, such as non-vortex state. Let us discuss 
here the local stability of our enumerated vortices. The 
stable ranges for the magnetization where the negative 
mode does not appear and system is stable for each vor- 
tex are listed below: 

(1) For the (1,0,-1) vortex configuration in the ferro- 



magnetic case, the stable range of M/N is determined 
as 

0.992 < M/N < 0.999 (ferromagnetic (1,0,-1) vortex) 

(34) 

The typical density profile and excitation spectrum are 
shown in Figs. 6(a) and 6(b) where the left hand side 
column depicts the density distributions for each com- 
ponent and the right hand side column shows the corre- 
sponding excitation spectrum as a function of the angu- 
lar momentum qg. The magnetization is M/N ~ 0.998 
here. 

(2) For the (1,0,-1) vortex in the antiferromagnctic 
case, the stable range is given by 

0.988 < M/N < 0.999 (antiferromagnetic (1,0,-1) vortex ) 

(35) 

We show the density profile and the excitation spectrum 
in Figs. 6(c) and 6(d). The magnetization is M/N = 
0.994 here. Comparing the above ferromagnetic case, 
because the magnetization decrease the 0-component fills 
in more the vortex core region. It pushes up the lowest 
excitation and stabilizes this vortex as a result. 

(3) In the Alice vortex of the antiferromagnetic case, 
the stable range is determined as 

0.972 < M/N < 0.999 (antiferromagnetic Alice vortex ) 

(36) 

We display the density profile and the excitation spec- 
trum for this stable Alice vortex in Figs. 6(e) and 6(f). 
The density distribution of the antiferromagnetic Alice 
vortex is very similar to that in the (1,0,-1) vortex shown 
in Fig. 6 (a) and (c). 

It is rather remarkable that these multiple components 
work together so as to stabilize the vortex locally. As 
the magnetization decreases from M/N = 1, the lowest 
eigenvalue at qg = — 1 becomes positive. Upon further 
decreasing beyond the critical magnetization, its eigen- 
value becomes negative again, indicating the local insta- 
bility. This was not seen in the single comp^nipt case 
where the vortex is always locally unstableoE3 This 
becomes stable only when the system is subjected un- 
der external rotation or when introducing the pinning 
potential. 

It is generally seen throughout these vortices that the 
component with the winding number fills in the core 
region of the +1 component with the winding 1. This 
0-winding component stabilize the vortex state. For the 
two component BECL-a similar stabilization mechanism 
is considered by Doi.LL§ 

The locally stable regions are limited to the cases 
where the magnetizations are close to the full polar- 
ized state. If the magnetization increases further to- 
wards the complete polarization, the system effectively 
becomes one component BEC and shows the lQfM insta- 
bility which was studied previously in details.BO This 
instability means the finite life time ^if, the vortex state 
which is observed experimentally! 
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Fig. 6. The density profiles (in the units of 10 20 m — 3 and fim); 
(a), (c) and (e), and the corresponding excitation spectra; (b), 
(d) and (f) for the locally stable vortices. The eigenvalues are 
normalized by hw and characterized by different marks, indicat- 
ing where the eigenvalues belong, (a) and (b) Ferromagnetic 
(1,0,-1) vortex, M/N = 0.998. (c) and (d) Antiferromagnetic 
(1,0,-1) vortex, M/N = 0.994. (e) and (f) Antiferromagnetic 
Alice vortex, M/N = 0.991. 



§6. Summary and Conclusion 

We have established a general framework to calculate 
various physical properties of a vortex in a cylindrical 
symmetric situation for a three component BEC sys- 
tem with rap — 1,0,-1 of F = 1. This allows us to 
evaluate the density distributions of the multiple compo- 
nent BEC, and the spatial variations of the triad (m, n, I) 
which fully characterize a vortex. This theoretical frame- 
work yields not only the above static properties, but also 
the dynamical properties, such as the collective modes. 
The latter is used to examine the local stability of a vor- 
tex. 

In the present paper, we have taken up three typical 
vortices which might be realized in the actual experimen- 
tal situations, namely the phase windings of the wave 
functions for the three components <f>i,<po and </>_i are 
given by (1,0,-1), (1,1,1) and (1, none, 0) as the repre- 
sentative vortices. 

We have characterized these vortices for both ferro- 
magnetic and antiferromagnetic interaction cases. In 
general, the 0- winding component, if any, fills in the core 



region. The densities of each component vary smoothly 
in the harmonically trapped confinement potential along 
the radial direction. The triad (m, n, 1) yields a par- 
ticular texture for each vortex. Especially the Z-vector 
which is proportional to the magnetization gives rise to 
a texture pattern, which should be observed by various 
experimental techniques. 

The Bogoliubov equations: eqs. (28) and (29) which 
yield the collective excitations have been extended to the 
three components. We calculate the excitations for the 
above three kinds of the vortices and examine the sta- 
bility condition as a function of the magnetization which 
can be controlled experimentally. Contrary to the sin- 
gle component BEC, the multi-component BEC allows 
a narrow intrinsic stability region for the magnetization, 
where the multiple components work to stabilize the vor- 
tex to raise the otherwise negative eigenvalue at qg = — 1 
to a positive value, that is, the non-winding component 
fills in the core region, effectively serving as a pinning 
potential. This self-organization in the multiple com- 
ponent systems for the intrinsic vortex stabilization is 
rather remarkable, because in the single component BEC 
the vortexas never stabilized without the external forced 
rotationBOJ) We naively expect that under the forced 
rotation these vortices enumerated here are easier to sta- 
bilize compare to the single BEC vortex. We will study 
it near future. 



[1] M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. Wieman 

and E. A. Cornell: Science 269 (1995) 198. 
[2] C. C. Bradley, C. A. Sackett, J. J. Tollett and R. G. Hulet: 

Phys. Rev. Lett. 75 (1995) 1687. 
[3] K. B. Davis, M.-O. Mewes, M. R. Andrews, N. J. van Druten, 

D. D. Durfee, D. M. Kurn and W. Ketterle: Phys. Rev. Lett. 

75 (1995) 3969. 

[4] M. R. Matthews, B. P. Anderson, P. C. Haljan, D. S. Hall, 
C. E. Wieman and E. A. Cornell: Phys. Rev. Lett. 83 (1999) 
2498. 

[5] B.P. Anderson, P.C. Haljan, C.E. Wieman and E.A.Cornell: 

Phys. Rev. Lett. 85 (2000) 2857. 
[6] K.W. Madison, F. Chevy, W. Wohlleben and J. Dalibard: 

Phys. Rev. Lett. 84 (2000) 806. 
[7] K.W. Madison, F Chevy, W. Wohlleben and J. Dalibard: 

|cond-mat/0004037 



[10] 

[11] 
[12] 
[13] 

[11] 

[15] 
[16] 



F. Chevy, K.W. Madison and J. Dalibard: Phys. Rev. Lett. 85 
(2000) 2223. 

T. Isoshima and K. Machida: J. Phys. Soc. Jpn. 66 (1997) 
3502. 

T. Isoshima and K. Machida: J. Phys. Soc. Jpn. 68 (1999) 
487. 

T. Isoshima and K. Machida: Phys. Rev. A 59 (1999) 2203. 

T. Isoshima and K. Machida: Phys. Rev. A 60 (1999) 3313. 

D.S. Rokhsar: Phys. Rev. Lett. 79 (1997) 2164. D. A. Butts 
and D. S. Rokhsar: Nature 397 (1999) 327. 

R. J. Dodd, K. Burnett, M. Edwards and C. W. Clark: Phys. 
Rev. A 56 (1997) 587. 

K. Doi and Y. Natsume: J. Phys. Soc. Jpn. (2001) 70 in press. 

J. Stenger, S. Inouye, D. M. Stamper-Kurn, H.-J. Micsncr, 
A. P. Chikkatur and W. Ketterle: Nature 369 (1998) 345. H.- 
J. Miesner, D. M. Stamper-Kurn, J. Stenger, S. Inouye, A. P. 
Chikkatur and W. Ketterle: Phys. Rev. Lett. 82 (1999) 2228. 
D. M. Stamper-Kurn, H.-J. Miesner, A. P. Chikkatur, S. In- 
ouye, J. Stenger and W. Ketterle: Phys. Rev. Lett. 83 (1999) 



8 



Author Name 



661. 

[17] For review of the spinor BEC. see D. M. Stamper-Kurn and 



W. Ketterle: cond-mat/0005001 



[18] T. Ohmi and K. Machida: J. Phys. Soc. Jpn. 67 (1998) 1822. 

[19] T.-L. Ho: Phys. Rev. Lett. 81 (1998) 742. 

[20] D. Vollhardt and P. Woffle: The Superfluid Phases of He 3 

(Taylor-Francis, London, 1990) 
[21] S.-K. Yip: Phys. Rev. Lett. 83 (1999) 4677. 
[22] U. Leonhardt and G. E. Volovik: JETP Lett. 72 (2000) 46. 



[23] H.T.C. Stoof: cond-mat/0002375 



[24] Karl-Peter Marzlin, Weiping Zhang and Barry C. Sanders: 

Phys. Rev. A 62 (2000) 013602 
[25] T. Isoshima, K. Machida and T. Ohmi: Phys. Rev. A 60 

(1999) 4857. 

[26] T. Isoshima, M. Nakahara, T. Ohmi and K. Machida: Phys. 

Rev. A 61 (2000) 063610. 
[27] M. Nakahara, T. Isoshima, K. Machida, S. Ogawa and T. 

Ohmi: Physica B 284-288 (2000) 17. 
[28] T. Isoshima, T. Ohmi and K. Machida: J. Phys. Soc. Jpn. 69 

(2000) No.12. 

[29] J. Martikainen, private communication. 



